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ABSTRACT 


Hueckel's  edge  detector  finds  the  best-fitting  ideal 
step  edge  to  a  given  picture  neighborhood,  by  expanding 
the  neighborhood  and  step  edge  in  terms  of  a  set  of 
nine  basis  functions.  The  simplest  case  of  this  ap¬ 
proach  uses  a  2-by-2  neighborhood  and  three  basis  func¬ 
tions.  This  case  is  solved  explicitly  using  elementary 
methods.  The  magnitude  of  the  best-fitting  step  edge 
for  the  neighborhood 


turns  out  to  be  the  Roberts  operator  max  ( | A— D | , ( B— C | ) 
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1.  Introduction 


Hueckel  [1-2]  developed  an  approach  to  edge  detection  based 
on  fitting  an  ideal  step  edge  to  a  given  picture  neighborhood. 
The  fitting  was  done  by  expanding  both  the  step  edge  and  the 
neighborhood  in  terms  of  a  set  of  orthogonal  basis  functions, 
and  minimizing  the  sum  of  the  squared  differences  between 
corresponding  coefficients.  To  simplify  the  computation,  the 
expansion  is  truncated;  Hueckel  used  a  nine- term  expansion. 

Several  simplifications  of  Hueckel 's  approach  have  also 
been  investigated.  Nevatia  [3]  used  a  subset  of  Hueckel' s 
basis;  O' Gorman  [4]  used  a  set  of  two-dimensional  Walsh  func¬ 
tions  defined  on  a  square;  Mero  and  Vassy  [5]  used  only  two 
basis  functions,  defined  by  diagonally  subdividing  a  square, 
to  determine  edge  orientation;  and  Hummel  [6]  used  a  set  of 
optimal  basis  functions  derived  from  the  Karhunen-Loeve  expan¬ 
sion  of  the  local  image  values. 

In  this  correspondence  we  present  an  elementary  treat¬ 
ment  of  edge  fitting  in  the  simplest  possible  nontrivial  case. 

AB 

We  use  a  2-by-2  picture  neighborhood  CD  ,  and  a  step  function 
s(x,y)  passing  through  the  center  of  this  neighborhood  (which 
we  take  to  be  the  origin) ,  defined  by 

fa  if  x  sin  8  >  y  cos  6 
s(x,y)  =  ^ 

^  b  otherwise 

We  use  only  three  basis  functions,  namely 


11  -1-1  1-1 

11  11  1-1 

It  turns  out  that  the  magnitude  of  the  best-fitting  step 
edge  derived  in  this  way  is  just  max  ( | A-D | , | B-C | ) ,  the  Roberts 
operator  [7].  Thus  this  correspondence  has  two  purposes:  to 
illustrate  how  Hueckel-type  edge  detectors  can  be  derived  using 
elementary  methods,  and  to  provide  a  new  motivation  for  the  max 
version  of  the  Roberts  operator. 
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2.  Derivation 

A  B 

Let  the  coefficients  of  f(x,y)  =  c  D  with  respect  to 
these  basis  functions  be  fg,f^,  and  f2»  respectively;  then 
readily  we  have 

fQ  =  (A+B+C+D) / 4  =  S/4 

f1  =  (-A-B+C+D)  /  4  (1) 

f2  =  (A-B+C-D) /4 

Similarly,  let  the  coefficients  of  s(x,y)  be  s0,Sj,  and 
s2?  then  readily 

Sq  =  (a+b)/2 

3l  ~  (b-a)  +  (a-b)  =  -  (b-a)  if  0  <  0  5  n 

=  (a-b)  +  (b-a)  *  3n2~"9'  (b-a)  if  it  <  0  <  2tt 

n  „  (2) 

s  '  =  —  (b-a)  +  —  (a-b)  *  £  (b-a)  if  -  l  <  0  <  y 
2  2n  2tt  "  2  2 

eJL  3n_0 

-  TT  (a"b)  +  TT  (b"a)  *  HT  (b"a>  if  7  £  e  -  T* 

2  2  2  2 
We  now  want  to  minimize  E  e  (fQ-s0)  +  (f^-s^  +  (f2~s2)  ; 

because  of  the*  way  and  s2  are  defined,  this  must  be  done 

separately  for  0  in  each  quadrant.  Actually,  since  a  and  b 

are  interchangeable,  by  symmetry  it  suffices  to  treat  the  first 

and  second  quadrants.  In  fact,  in  (2  ),  if  we  replace  0  by 

0+tt  and  interchange  a  and  b  in  3w^29  (b-a),  we  obtain  (b-a) 

and  similarly  2— ’  (b-a)  yields  —  (b-a) . 

n  it 


Also,  substituting  S/4  *  (a+b)/2  and  (  5 )  in  either  equation 
of  ( 4  )  gives 


(b-a)  =  0 


(7) 


-A-B+C+D  20-rr 

3  7tT 

so  that  we  also  have 


A-B+C-D 

- 3 - 


§  (»-a) 


0 


(8) 


If  9  is  not  0  or  n/2,  and  b  ^  a,  we  can  divide  (  7  )  by 
(  8 )  (or  vice  versa)  to  obtain 


A-B+C-D  _  9 

-A-B+C+D  0  rr  (9) 

e"2 

from  which  we  readily  have 


0 (2A-2D)  =  j  (A-B+C-D) 


or  (if  A/D) 


_  £  t  A-B+C-D  _  n  .  ,  _  B-C  , 
“  4  A-D  “  7  ^ 1  "  A-D 


(10) 


Combining  this  with  (  8  )  gives  b-a  =  A-D,  and  combining  this 
with  b+a  =  S/2  gives 


a  -  (-A+B+C+3D) / 4  «  S/4  -  (A-D)/2 
b  =  ( 3A+B+C-D) / 4  =  S/4  +  (A-D)/2 


(ID 


Note  that  by  (71,  (8  ),  and  the  fact  that  a+b  *  S/2,  we  actually 
2 

have  E  «  0  for  this  solution.  It  is  easily  verified  that  if 
we  assume  6  »  0  or  0  *  ^  in  (3),  and  set  the  partial  deriva¬ 
tives  with  respect  to  a  and  b  equal  to  zero,  we  obtain 


special  cases  of  this  solution.  In  fact,  for  0  =  0,  we  find 
that  A+C  »  B+D,  b  =  (A+B)/2,  and  a  «  (C+D)/2;  while  for  0  =  7 
we  get  A+B  =  C+D,  b  =  (A+C)/2,  and  A  =  (B+D)/2. 

In  the  second  quadrant,  analogously,  we  get 

a  =  (A+3B-C+D) / 4  =  S/4  -  (C-B)/2 

(12) 

b  =  (A-B+3C+D) /4  =  S/4  +  (C-B)/2 


a  _  n  r ,  .  A-D  , 

6  =  4  1 3  +  B^C  ]  (13) 

It  can  be  verified  that  the  first  and  second  quadrant  solutions 
agree  if  6  =  tt/2.  Moreover,  note  that  for  (10)  to  actually  lie 
in  the  first  quadrant  we  must  have 

-1  <  *z£  <  1 
x  "  A-D  - 

which  is  evidently  equivalent  to  | B— C |  5  |A-D|.  Similarly, 
for  (13)  to  actually  lie  in  the  second  quadrant  we  must 
have 

-i  <  Sr£  <  1 

B-C 


which  is  evidently  equivalent  to  j  A— D  j  <  |B-C|.  Thus  by  compar¬ 
ing  the  magnitudes  and  signs  of  A-D  and  B-C  we  can  choose 
the  appropriate  best- fitting  step  edge  for  the  given  neighbor¬ 
hood  Note  that  if  A-D  =  B-C  we  have  0  =  0  in  (10)  and 

0  an  in  (13);  moreover,  in  this  case  (11)  and  (12)  also  agree, 
with  a  and  b  interchanged.  Similarly,  if  A-D  =  C-B,  we 
have  0  =  |  in  both  (10)  and  (13),  and  here  (11)  and  (12)  agree 


too 


In  summary,  the  "best-fitting"  step  edge  to  “  is  found 
as  follows: 

If  |  B-C  |  <  I A-D  |  ,  then  ©  =  J  [  1  -  f=§  ]>  and  a,b  are 

given  by  (11) 


If  |  B-C  |  >  |  A-D  | ,  then  6  =  j  [  3  +  ],  and  a,b  are 

given  by  (12) 

The  magnitude  |a-b|  of  the  edge  is  | A— D |  in  the  first  case, 

and  | B-C |  in  the  second  case;  in  other  words,  the  magnitude 

is  max  ( | A-D | , | B-C | ) .  Note  that  this  is  just  the  magnitude  of 

the  Roberts  operator,  using  the  max  of  the  absolute  differences 

rather  than  the  square  root  of  the  sum  of  the  squares  [7]. 

(The  slope  0,  on  the  other  hand,  is  not  the  arc  tangent  of  the 

ratio  of  these  differences;  but  its  value  is  reasonable,  e.g., 

, -  AB  12 


IT 


3.  Conclusion 


We  have  presented  an  elementary  derivation  of  step  edge 
fitting  in  the  simplest  nontrivial  case:  a  2-by-2  neighbor¬ 
hood  and  three  basis  functions.  It  turns  out  that  the  magni- 
tude  of  the  best-fitting  edge  to  CD  is  max  ( | A-D | , | B-C | ) , 
which  is  a  commonly  used  version  of  the  Roberts  edge  detector; 
thus  our  derivation  provides  a  new  motivation  for  that  detector. 
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